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1  Summary 


The  inverse  scattering  solution  of  shape  and/or  material  parameter  reconstruction  is  of¬ 
ten  posed  as  a  problem  in  nonlinear  minimization  of  an  objective  function  with  respect 
to  N  (usually  large)  number  of  unknown  model  parameters  characterizing  the  scatterer. 
The  minimization  procedures  are  usually  iterative,  and  require  the  gradient  of  the  ob¬ 
jective  function  in  the  unknown  model  parameter  vector  in  each  stage  of  iteration.  For 
large  N,  finite-differencing  becomes  numerically  intensive,  and  an  efficient  alternative  is 
domain  differentiation  in  which  the  full  gradient  is  obtained  by  solving  a  single  scattering 
problem  of  an  auxiliary  field  using  the  same  scattering  operator  as  that  of  the  forward 
solution.  A  well-known  technique  in  this  direction  is  the  so-called  adjoint  field  method 
which  obtains  the  gradient  by  variationally  minimizing  an  augmented  objective  function 
that  includes  the  reduced  wave  equation  via  a  Lagrange  multiplier.  Results  are  reported 
mostly  for  compact  objects.  This  report  presents  the  domain  derivative  calculation  of  the 
gradient  for  a  one-dimensional,  locally  perturbed,  infinitely  long  dielectric  interface.  The 
method  is  non-variational,  and  algebraic  in  nature  in  that  it  evaluates  the  gradient  by 
directly  domain  differentiating  the  scattering  equations.  The  computations  are  straight¬ 
forward,  and  easy  to  follow.  The  mathematical  transformation  of  the  scattering  problem 
into  the  corresponding  problem  for  the  differentiated  fields  can  be  visualized  explicitly. 
The  formulation  of  and  the  motivation  behind  introducing  the  auxiliary  field  are  explic¬ 
itly  demonstrated.  Closed-form  analytic  expressions  are  obtained  for  the  gradients  for 
electromagnetic  transverse  electric  /  transverse  magnetic  (TE/TM)  scattering  from  di¬ 
electric  rough  surfaces,  and  for  scalar  wave  scattering  from  Neumann  and  Dirichlet  rough 
surfaces.  Results  compared  with  those  for  compact  scatterers.  Finally,  the  relationship 
between  our  results  and  Lorentz  reciprocity  is  pointed  out  and  clarified. 
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2  Introduction 


Consider  the  scattering  of  an  electromagnetic  field  from  a  one-dimensional,  infinitely 
long,  locally  perturbed,  flat  interface,  T,  shown  in  Figure  1.  It  is  assumed  that  the 
interface  is  located  at  x 2  =  0,  and  is  invariant  along  £3.  The  graph  of  the  perturbation, 
7  =  {xi,  £2  =  /(£  1)},  —a  <  £1  <  a,  is  assumed  to  be  smooth.  Moreover,  the  profile,  f{x  1) 
is  a  single- valued  function  of  X\.  The  surface  is  thus  not  re-entrant.  T  separates  R 2  into  an 
upper  half-plane,  h2+  =  {£i,£2  >  /(£i)},  and  a  lower  half-plane,  =  {£i,£2  <  f(x  1)}, 
with  wavenumbers,  k+ ,  k~ ,  and  dielectric  constants,  e+  and  e~ ,  respectively.  It  is  further 
assumed  that  the  perturbed  surface  is  illuminated  at  a  non- grazing  angle  ( \6inc\  <  vr/2) 
by  a  beam  (-0mc)  of  finite  width,  e.g.,  a  beam  with  a  Gaussian  angular  spectrum  [40]  or  a 
tapered  Thorsos  beam  [44]  having  negligible  energy  beyond  |£i|  >  A,  with  4>a.  Finite 
width  beams  and  their  properties  are  discussed  in  [23,  45,  24], 


Figure  1:  A  Schematic  Illustration  of  the  Scattering  Geometry 


The  following  equations  [45,  37,  34,  32,  8,  38,  41,  28]  describe  the  scattering  problem. 


(£+^+)(x)  = 
= 

i\)+  = 

= 

In  Eqs.  (l)-(4),  CA  =  A  +  (A:±)2,  A  being  tl: 
is  assumed.  ^  are  the  total  fields,  and  ipz. 
The  unit  normal,  h.  is  given  by 


0,  x  G  fl+, 

(1) 

0,  x  G  G~, 

(2) 

V’-,  x  g  r, 

(3) 

a+jp-,  x  e  r. 

(4) 

:  Laplacian,  and 

a  time  harmonicity  of  e~lujt 

their  normal  derivatives  on  the  boundary. 


n  = 


y/1  +  (/'(*  1))2 


[f\x  1),  l]5 
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and  points  into  fi+.  f  is  the  derivative  of  /  in  x\,  and  the  superscript  T  denotes  trans¬ 
pose.  Moreover,  =  ^mc  +  ij)sc  is  the  sum  of  the  incident  (ipinc)  and  scattered  (?/Pc) 
held.  For  the  scattering  geometry  in  Figure  1,  ijjsc  satisfies  [34,  10,  9,  2]  the  radiation 
condition  at  infinity  [7].  The  parameters,  a1*1,  depend  on  the  material  constants  in  H1*1, 
and  the  scattering  process.  Applications  to  electromagnetic  TE  and  TM  scattering  from 
a  dielectric  surface,  and  for  scalar  wave  scattering  from  Neumann  and  Dirichlet  surfaces, 
are  considered  later. 

For  an  arbitrary  surface,  a  closed  form  solution  of  (l)-(4)  is  not  expected.  Several 
approximate  asymptotic  techniques  have  been  developed  during  the  past  few  decades  for 
this  kind  of  problems.  Two  classical  approaches  are  the  small  perturbation  method  (SPM) 
[39,  3,  46],  and  the  Kirchhoff  approximation  (KA)  [4,  3,  46].  SPM  may  be  used  when 
the  surface  fluctuations  are  small  and  smooth.  KA  is  applicable  for  problems  where  the 
radii  of  fluctuations  are  large  compared  to  the  signal  wavelength.  One  observes  that  SPM 
leads  to  a  linear  relationship  between  the  scattered  held  and  the  function  describing  the 
surface  fluctuation.  Hence  one  can  readily  employ  the  Fourier  transform  technique  to 
invert  the  scattered  held  to  reconstruct  the  rough  interface  [48],  [1].  It  turns  out  that 
even  when  using  KA,  the  scattered  helds  may  be  manipulated  such  that  the  the  rough 
interface  can  be  reconstructed  using  the  Fourier  transform  technique  [47],  [42],  [30].  These 
classical  approaches  to  rough  surface  scattering,  viz.,  SPM  and  KA,  are  single  scattering 
approximations  and  hence  a  linear  relationship  between  the  scattered  held  and  the  under¬ 
lying  rough  surface  is  possible.  However  if  the  surface  roughness  is  not  small  and  smooth, 
multiple  scattering  becomes  important.  Such  processes  are  nonlinear  in  character  and 
hence  the  solution  to  (l)-(4)  is  formally  expressed  as:  i/;sc  =  F(f),  where  T  is  a  nonlinear 
operator  that  transforms  the  boundary,  T,  into  scattered  helds.  The  determination  of 
,ijjsc,  given  T  and  /,  constitutes  the  direct  or  forward  problem.  The  inverse  problem  is 
to  reconstruct  the  surface  prohle,  /,  where  T  and  i\)sc  are  known.  Formally,  the  inverse 
solution  is:  /  =  JF_1{^SC},  the  result  of  applying  the  inverse  scattering  operator,  to 
the  data,  {^sc}.  However,  because  of  the  nonlinear  nature  of  T ,  the  inverse  problem  is 
frequently  solved  [7,  15,  16]  by  transforming  it  into  a  problem  in  nonlinear  optimization. 
The  solution  is  obtained  as  the  minimizer  of  an  objective  function,  77(/),  defined  as: 

1  M 

3(5 )  :=  5 

m=  1 

M 

=  5^(Ay)-(fm)(AiV.)(fm).  (5) 

771—1 

In  (5),  A 'V(Tm)  :=  —  i/jmeas(xm)  is  the  data  misfit  function,  i.e.,  the  difference 

between  the  measured  meas )  and  a  theoretically  predicted  scattered  held  (/ippr),  both  at 
a  detector  location,  xm.  The  complex  conjugation  is  denoted  by  *,  and  M  is  the  total 
number  of  the  detectors.  The  superscript,  5,  on  A5  indicates  that  meas  is  noisy,  5  being 
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the  level  of  the  measurement  noise.  Furthermore,  exact  as  5  — >  0  [13,  27],  where 

ijjs  is  a  regularized  solution. 

Let  p  =  {pi}i  be  the  parametrization  of  7.  The  inverse  solution  is  then  given  by: 
p  =  argminXj  J (p) ,  Xp  being  the  space  of  admissible  parameters.  The  inverse  problem 
being  ill-posed,  regularization  [13]  is  necessary  in  order  to  obtain  a  meaningful  solution. 
The  standard  procedures  for  the  minimization  of  (5)  are  iterative  in  nature,  and  require 
that  the  gradient  of  the  objective  function  (or  the  scattered  fields)  in  the  unknown  model 
parameter  vector  be  obtained  in  each  stage  of  iteration.  Let  this  derivative  of  J(f)  be 
denoted  by  dj{f). Then: 

M 

dj{f)  =  -Re^(A  S,4>)*{xm)'4>s\xm).  (6) 

771=1 

ips'  in  (6)  is  the  domain  derivative  [26,  18,  19],  of  fi>sc,  that  is,  the  derivative  of  the  scattered 
held  with  respect  to  the  variation  in  /.  Equation  (6)  can  also  be  written  as: 

dJ(f)  =  -RelJ^A'V)*},  (7) 

where  JT  is  the  transpose  of  the  Jacobian,  J. 

The  procedure  most  often  employed  in  inverse  problems  for  the  computation  of  deriva¬ 
tives  of  the  cost  function  w.r.t.  parameters  of  the  problem  is  numerical  [17].  For  instance, 
in  our  problem,  one  can  use  a  N  parameter  model  to  represent  the  rough  interface  7 
and  employ  finite  difference  method  to  compute  the  domain  derivative.  This  is  a  very 
general  method  applicable  any  complex  problem.  However,  differentiating  fi>sc  with  re¬ 
spect  to  a  parameter  vector  of  length  N  involves  solving  N  distinct  scattering  prob¬ 
lems  for  N  distinct  scatterers.  If  N  is  large  (which  is  often  the  case  in  practice),  then 
clearly,  the  finite-difference  method  becomes  numerically  intensive,  besides  accumulating 
large  round-off  errors.  It  essentially  involves  the  calculation  of  the  large  Jacobian  ma¬ 
trix,  J.  However,  (7)  shows  that  the  evaluation  of  the  derivative  involves  evaluating  the 
product  term,  Jr(A5^)*,  and  not  the  Jacobian  explicitly.  In  the  adjoint  field  method 
[31,  12,  11,  35,  36,  43,  22,  14],  this  product  is  obtained  as  the  solution  of  the  forward 
problem  in  Eqs.  (l)-(4),  but  with  a  source  term  as  in  the  R.  H.  S.  of  (6).  The  source 
is  obtained  by  backpropagating  the  data  misfit  functions  to  the  scatterer  space.  The 
boundary  value  problem  for  the  adjoint  held  is  derived  by  applying  variational  techniques 
to  minimize  an  augmented  objective  function,  namely,  J  in  (5)  plus  a  constraint  term  in¬ 
volving  the  reduced  wave  equation  via  Lagrange  multiplier.  The  upshot  of  all  this  is  that 
the  entire  gradient  is  evaluated  by  solving  only  one  single  scattering  problem  in  which  the 
same  operator  as  that  of  the  forward  scattering  problem  is  used.  An  important  feature  of 
this  method  is  that  it  is  essentially  independent  of  N.  There  are  two  different  ways  one 
can  carry  out  the  domain  derivative  calculation  using  the  adjoint  method.  Depending  on 
the  numerical  procedure  used  for  solving  the  forward  problem  one  can  use  appropriate 
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discretization  and  carry  out  the  numerical  calculation  of  the  derivatives.  This  is  a  rela¬ 
tively  more  general  procedure.  In  some  problems  it  may  be  possible  to  obtain  analytical 
closed-form  expression  for  the  domain  derivative.  This  not  only  offers  valuable  physical 
insight  into  the  problem  but  also  gives  the  flexibility  in  the  choice  of  algorithm  used  for 
shape  reconstruction. 

In  this  report,  we  focus  attention  on  obtaining  closed-form  expression  for  the  domain 
derivatives.  There  are  several  distinguishing  features  in  the  calculations  of  this  report.  In 
the  literature  we  notice  that  domain  derivatives  of  the  scattered  holds  have  been  calcu¬ 
lated  mostly  for  compact  scatterers,  i.e.,  scatterers  occupying  a  finite  volume  of  space.  In 
this  report,  on  the  other  hand,  an  infinitely  long,  locally  perturbed,  dielectric  interface  is 
considered.  The  boundary  variation  in  this  case  is  relatively  more  subtle  than  that  for  a 
compact,  volumetric  object.  Notice  that  we  use  neither  variational  procedures  nor  an  aug¬ 
mented  objective  function  in  this  report.  The  calculations  are  purely  algebraic  in  nature 
in  that  the  scattering  equations  are  directly  domain-differentiated  without  any  interven¬ 
ing  mathematical  formalism.  In  this  manner  of  presentation,  the  calculational  procedures 
are  found  to  be  straightforward,  physically  transparent,  and  easy  to  follow.  Moreover,  the 
details  of  the  mathematical  transformation  of  the  scattering  problem  into  the  correspond¬ 
ing  problem  for  the  differentiated  fields  can  be  traced  explicitly.  The  present  calculations 
also  invoke  an  auxiliary  scattering  problem,  as  in  the  adjoint  held  case.  However,  the 
mathematical  formulation  of  and  the  motivation  behind  the  auxiliary  problem  are  more 
transparent  here  than  in  the  adjoint  held  calculations  where  the  problem  is  a  result  of  the 
mathematical  manipulations.  Also,  the  gradients  are  calculated  for  the  electromagnetic 
TE  and  TM  scattering  from  a  dielectric  rough  surface  as  also  for  scalar  wave  scattering 
from  a  Neumann  and  a  Dirichlet  surface.  Finally,  the  relation  of  the  results  thus  obtained 
with  the  Lorentz  reciprocity  is  discussed  and  clarified. 

The  report  is  organized  as  follows.  Some  important  relations  from  domain  differential 
calculus,  relevant  to  the  present  calculations,  are  summarized  in  Section  3.  These  are  then 
used  in  Section  4  to  derive  the  boundary  value  problem  for  the  derivative  helds,  ■0±/,  of 
riJj±.  The  differential,  dj (/),  and  the  gradient,  V/Jr(/),  of  ),  are  evaluated  in  Section 
5.  The  method  is  next  applied  to  electromagnetic  TE  and  TM  scattering,  and  Neumann 
and  Dirichlet  boundary  conditions  in  Section  6.  Discussions  are  presented  in  Section  7, 
followed  by  a  summary  in  Section  8. 
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3  Some  Relevant  Identities  From  Domain  Differentiation 

Let  U  be  a  piecewise  smooth  bounded  open  region  in  R 2  with  boundary  E.  A  continuously 
differentiable  vector  field,  V,  is  defined  in  U.  It  is  assumed  that  V  or  its  normal  derivative 
vanishes  on  £.  Consider  the  mapping,  Tr(V(x))  =  I  +  tV(x)  +0(t2),  feW,  continuously 
mapping  U  onto  itself.  /  is  the  identity  operator,  and  the  parameter,  r  G  [0,2],  2  G  R1.  Tr 
introduces  a  change  of  variable  from  x  — ►  xT  =  x  +  rV(x)  +  0(r2).  Let  12  C  U  be  a 
domain  in  U.  Then  Tr(Q)  =  I2r,  where  f2r  =  {xT\xT  =  x  +  rV(x),  VT  G  f2}.  Moreover, 
C  W,  Vr  G  [0,  t].  r  defines  the  magnitude  of  the  deformation.  In  fluid  dynamical 
problems,  r  plays  the  role  of  the  time  variable.  These  are  illustrated  in  Figure  2.  It  can 
be  shown  [49]  that  for  small  r,  Tt  is  a  diffeomorphism  in  U ,  that  is,  a  bijection  with  a 
continuous  inverse,  with  same  degree  of  smoothness  as  the  vector  field.  This  inverse  is 
given  by:  T~\V)  =  I  -  rV  +  0(r2). 


Figure  2:  Domain  Variation  for  a  Compactly  Supported  Scatterer 

Let  (f){x)  be  the  solution  of  a  boundary  value  problem  defined  on  12,  which  is  considered 
to  be  the  unperturbed  domain,  i.e.,  12  =  f2r=o.  Let  Tr  transform  12  to  the  perturbed 
domain,  Tr(i 2)  =  I2r,  and  let  <j)T(xT),  xT  G  I2r,  be  the  corresponding  solution,  now  defined 
on  I2t.  Note  that  r  appears  twice  in  (j>T(xT),  once  as  a  subscript  in  </>T,  and  second  time  in 
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xT,  which  is  a  running  variable.  It  is  the  dependence  of  (p  on  r  that  gives  rise  to  domain 
differentiation  in  place  of  the  regular,  partial  derivatives.  The  Eulerian  (also  called  total 
or  material )  derivative  of  <p  at  x,  denoted  by  <p(x),  is  defined  as: 

4>{x)  =  lim  —  [cj)T{xT)  —  4>(x)\.  (8) 

r^O  T 

The  derivative  in  (8)  can  also  be  written  as  the  directional  derivative: 

4>{X)  =  ^(0  •  Tr)|r=0 

dr  (t" j  3't)  | r=0 •  (9) 

<j),  as  defined  in  (8)  (or  (9)),  is  a  Gateaux  differential  [29],  which  is  defined  in  a  pointwise 
sense  at  x  e  fb  Equation  (9)  can  be  further  expressed  as: 

<j)(x)  =  lim  ~[<i>T(xT)  -  <f>(x)\, 

t  — >0  T 

=  lim  —[(j)T{x)  —  <f>(x)]  (10) 

r^O  T 

+  lim  -  [(j)T(xT)  -  4>t(x)\. 

t—> 0  T 

The  solution  is  assumed  to  be  defined  in  U  x  /  C  Rn+1,I  =  [0,  t]  C  Rl .  0r  is  then 
considered  to  be  the  restriction  of  the  solution  defined  on  flT.  Provided  that  the  limit  on 
the  first  term  in  the  R.  H.  S.  in  (10)  exists,  the  equation  can  be  recast  as: 

<fi(x)  =  4>'(x)  +  V(x)  ■  V<p(x),  (11) 


where 

4>'{x)  =  lim  —[(j)T{x)  —  <t>(x)\  (12) 

r^O  T 

is  the  partial  derivative  of  <t>  at  x,  i.e.,  the  derivative  when  x  is  held  fixed.  The  term, 
V  ■  V0,  is  the  convective  part.  Furthermore,  since  r  and  x  are  independent  variables,  dT 
and  dx  commute,  yielding  the  following  important  relation: 

(V0)'  =  V<//,  (13) 


The  Eulerian  derivative  of  a  domain  functional  can  be  obtained  from  (12).  Let 

F((pT;  flT)  =  /  cpT(xT)  dxT 
J  nT 

be  a  domain  functional  of  0r  over  the  perturbed  domain,  h2r.  It  can  be  written  as: 


F((j)T]Q )  =  /  (j)T(x  +  rV (x))  \  Jt\  dx. 

Jn 


(14) 
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\Jt\  is  the  Jacobian  matrix  of  Tr,  and  ( Jr)ij  =  $ij  +  Vi,j,  Vi,j  —  dVi/dxj.  Direct  calculation 
yields:  \Jt\  =  V  ■  V.  Moreover,  Jt\t=o  =  I-  Then  differentiating  Eq.  (14)  in  r  yields  the 
Eulerian  derivative,  dF,  namely 

dF(f>]  Q)  =  /  <f\x)  dx  +  /  <j>(x)v(x)  ds.  (15) 

Jn  Jr 

v  —  V  ■  n  is  the  normal  component  of  V  on  T.  It  is  demonstrated  (Section  5)  that  dF(q f  Q) 
is  linear  in  the  perturbation.  The  Gateaux  derivative  [29],  therefore,  exists,  and  F((f>]Q) 
has  a  gradient. 

Let  us  further  introduce  the  following  identities. 


V0  =  Vr0  +  n(n  ■  V0)  =  Vr0  +  n0jTl, 

Vr0  •  Vpu  =  V</>  •  Vu  -  0>nun, 

Vr  •  (zVr(p)u!  ds  =  —  [  z(Vr4>  •  Vrw)  ds. 


(16a) 

(16b) 

(17) 


Vr  =  V  —  h<9n  is  the  surface  gradient,  and  the  functions  w  and  z  are  defined  in  a 
neighborhood  of  T.  Eq.  (17)  is  sometime  referred  to  as  Stokes'  identity  which  can  be 
established  using  differential  geometry  [33].  Suffice  it  to  say  that  the  surface  gradient  of 
a  function  u  defined  on  the  boundary  is  the  restriction  to  the  boundary  of  the  lifting  of 
u  defined  in  the  tubular  covering  of  the  boundary.  Finally,  we  introduce  Green’s  second 
theorem  [25,  7],  which  is: 

/  +  V0  •  V"0]  dx  —  f'fin  ds.  (18) 

Jn  Jr 

The  unit  normal,  n,  is  directed  in  the  exterior  of  the  domain  fl.  <fT{xT )  is  a  domain 
functional. 

The  above  discussions  were  for  a  closed,  bounded  domain  (Figure  2).  In  the  case  of 
the  locally  perturbed  interface  in  Figure  1,  the  deformation  is  modified  to  that  shown  in 
Figure  3.  As  earlier,  the  velocity  held  is  assumed  to  vanish  on  E.  The  surface  profile, 
f(x i),  separates  U  into  U+  and  U~ .  Thus  U  =  U+  +U~.  As  the  curve  f(xi)  is  deformed 
into  fT(x i),  U±  are  deformed  to  Uf,  but  the  overall  domain,  U ,  remains  unchanged.  The 
relations  which  were  derived  above,  must  now  be  applied  to  the  regions,  Uf  and  Uf .  All 
derivations  assume  C'2-smoothness  of  the  boundary  function.  Thus  6  C2(fl  ),  and 
the  deformation  operator,  Tr  is  C'2-diffeomorphic.  Furthermore,  it  is  known  that  for  the 
smoothness  condition  assumed  here,  the  solutions  can  be  extended  to  the  neighborhood  of 
the  boundaries  of  D,  and  VT.  For  details,  we  refer  to  [20].  We  also  note  that  the  application 
of  Green’s  theorem  to  U  in  Figure  3  may  be  problematic.  However,  the  volume  integrals 
involving  A-f/A,  the  Laplacian  appears  as  being  multiplied  by  functions  that  are  smooth 
in  the  domain,  and  vanishes  on  the  boundary.  The  application  of  Green’s  theorem  then 
does  not  cause  problems.  Next  we  derive  the  BVP  for  the  domain  derivatives  of  in 
Eqs.  (l)-(4)  using  the  relations  that  were  presented  in  this  section. 
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4  The  scattering  problem  for  the  field  derivatives 


In  this  Section,  we  determine  how  the  scattering  equations  (l)-(4)  for  the  unperturbed 
surface  are  transformed  as  the  surface  is  deformed  by  Tt.  Scattering  in  the  transformed 
domains,  U± ,  is  described  by: 


(Cyv+)(£r) 

=  0, 

XT  G 

(19) 

(Cyv)(£r) 

=  0, 

xT  eWT_, 

(20) 

= 

Vv,  Xt 

£  7 T1 

(21) 

a+C„, 

xT  e  7 T. 

(22) 

The  Helmholtz  equations  (19)  and  (20)  are  transformed  first.  The  transformation  of  the 
boundary  conditions  in  Eqs.  (21)  and  (22)  follows  next.  From  the  discussion  at  the  end 
of  the  previous  Section,  Tr  =  yT,  and  all  differentiations  are,  therefore,  with  respect  to  7. 
Taking  the  derivatives  of  Eqs.  (19)  and  (20),  and  using  (11)  yields: 

(A ^Y  +  V  -V(A ^±)  +  {k±)\(^±)' +  V  ■V(r/;±))  =  0. 

Now,  (A^)'  =  A (t^y  by  Eq.  (13),  and  £±^±  =  0  in  W±  by  virtue  of  Eqs.  (1)  and  (2). 
It  then  follows  that: 


(£+r/A')(f)  =  0,  xeU+,  (23) 

(/ryr')(£)  =  0,  xeu~.  (24) 

Now,  in  the  exterior  of  U,  =  '0=b/,  since  the  observation  point  is  held  fixed.  The 
convective  term  then  vanishes,  and  C±'ij>±'  vanishes  outside  IA.  Therefore, 

C+ip+\x)  =  0,  (25) 

C~rp~'(x)  =  0,  x  G  .  (26) 


Exactly  similar  considerations  are  applied  to  the  transformation  of  the  first  boundary 
condition  in  Eq.  (21).  The  result  is: 

yH'-yr'  =  -f  •  v(yA  -  yo 

=  -V  ■  v7(v>+  -  r )  -  »W+  - 1*,;). 

=  (a  -  1)V>>.  (27) 


In  Eq.  (27),  V7  is  the  tangential  derivative  on  7.  Furthermore,  a  =  a~/a+,  and  the 
boundary  conditions  in  Eqs.  (3)  and  (4)  were  used  in  arriving  at  the  result. 

The  Eulerian  derivative  of  the  second  boundary  condition  in  Eq.  (4)  is  calculated  as 
follows.  First  we  introduce  the  following  notations. 

r±  r  r 


$±dx=  $+dx±  dx 


'v± 


'v+ 


IV- 
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Figure  3:  Domain  Variation  for  the  Locally  Perturbed  Surface  of  Figure  1 


for  the  volume,  and 


J  d>±  ds  =  J  ($+  ±  $")  ds 


for  the  surface  integral.  Also,  define  a  regular  function,  u  G  hi.  which  vanishes  on  £,  and 
consider  the  domain  functional: 


r-  + 


F(U+,Ur  )  =  /  (aT dxT. 


Apply  Green’s  theorem  (Eq.  (18))  to  Eq.  (28),  and  obtain: 


(28) 


F(K,l C)  = 


{aT^n}  u  ds7 


—  /  Vr  )  •  VTw  dxT  +  /  dfT.  (29) 

Ju±  JuF 

Vr  is  the  gradient  with  respect  to  xT,  and  the  negative  sign  on  the  qT  integral  is  due  to 
the  opposing  signs  of  the  n  on  the  surface.  The  L.H.S.  in  Eq.  (29)  vanishes  because  of 
Eqs.  (19)  and  (20).  Moreover,  since  oTi\)+n  =  ar^~n  from  the  boundary  condition  in  Eq. 
(22),  the  surface  integral  in  Eq.  (29)  vanishes  yielding: 


-+ 


tu? 


Vr  (a=F,0T  )  '  ^tU  dxT  +  /  {^(A;111)2}  ip±u  dxT  =  0. 


(30) 


'W 


We  now  take  the  Eulerian  derivative  in  Eq.  (30),  use  Eq.  (11),  and  simplify.  The  result 


is: 


/w± 


[V(a=F,0±/)  •  Vu  —  a¥{k±)2'^±'u\  dx  = 
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[V(a=F'^=t)  •  Vu  —  aT (k±)2'ijju  v  ds. 


(31) 


We  used  the  fact  that  7+  =  7”  =  7,  and  ^ +  =  —  tp  on  7  by  the  boundary  condition 

in  Eq.  (3).  We  now  replace  (AA)2^1*1'  by  —A^'  (cf.  Eqs.  (23)  and  (24))  in  Eq.  (31)  to 
obtain: 

+ 

[A (a^'V  +  V(aT'0±/)  -Vu]  df  = 


L 


—  /  [V(a=F'^±)  ■  Vtt  —  |aT(A;±)2  j  v  ds.  (32) 

J  7 

Applying  Green’s  theorem  to  the  L.H.S.  in  Eq.  (32)  yields: 

f  [A(a:=F'0±/)'u  +  V(a:=F'0±,)  ■  V«]  dx  —  f  (o'=F^^/)  u  ds.  (33) 

J  IA^~  7 

The  Erst  surface  integral  in  the  R.H.S.  in  Eq.  (32)  can  be  evaluated  as  follows. 

/  [V(a=F'^±)  •  Vu]  v  ds  —  /  [V7(a=F'0±)  ■  V7m  +  «,„]  r  ds.  (34) 

J  7  J  7 

As  before,  ct-^  =  o.+^~n  from  the  boundary  condition  in  Eq.  (4).  The  second  integral 
in  the  R.H.S.  in  Eq.  (34),  therefore,  vanishes,  yielding: 

/  [V(aT^)  ■Vu]vds  = 


[V7(a=F'^d=)  •  V7w]  v  ds 

~)  f  {(V7^+)  •  (V7u)}v  ds 
=  (a+  —  a~)  f  {V7  •  (uV7'0~l~)}  u  ds, 


=  [a  —a 


(35) 


in  which  the  identity  in  Eq.  (17)  was  used.  Moreover,  V7^  =  V7^+.  Replace  Eqs.  (33) 
and  (35)  in  Eq.  (32).  This  results  in: 


'7 


{afFV’,re/}  u  ds  =  (a  —  a+)  /  V7  •  (uV7'0+)u  ds+  [a:=F(A;±)2'0u]  u  ds, 

t/7  J  7 

from  which  it  follows  that: 

oi~ ~  a+^n  —  («_  —  ct+)  {V7  •  (uV7^+)}  +  (a_A;+2  —  a+k^2)^v.  (36) 

Equation  (36)  establishes  the  boundary  conditions  on  7  for  the  normal  derivatives,  t/W, 
of  the  derivative  fields, 
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Collecting  the  results  in  Eqs.  (25),  (26),  (27)  and  (36)  together  yields: 


C+^+\x) 

=  0,  x  g  n+, 

(37) 

=  0,  x  G  fl_, 

(38) 

i/j+l  - 1\)~' 

=  (a-  l)ip%v,  x  G  7, 

(39) 

i+/  -  a+'ib~/ 

n  ^  y  ,n 

=  (cT  —  ct+)V7  •  (uV7'0+) 

+(a~k+ 2  —  a+k~2)^+v,  x  G  7, 

(40) 

^+/ 

=  x  e  r\7, 

(41) 

=  a+ip^,  x  e  r\7. 

(42) 

The  last  two  boundary  conditions  in  Eqs.  (41)  and  (42)  follow  simply  from  the  fact  that  on 
the  flat  part  of  T,  the  observation  points  are  held  constant.  Equations  (37)  through  (42) 
constitute  the  scattering  problem  for  the  Eulerian  derivatives  of  the  Helmholtz  problem 
in  Eqs.  (l)-(4).  As  can  be  seen  from  the  above  equations,  the  scattering  problem  for  the 
derivative  fields  is  essentially  the  same  as  that  for  the  unperturbed  surface  except  that 
the  boundary  conditions  on  7  for  '0±/  now  depend  upon  the  deformation  held  through  the 
parameter,  v.  The  deformation,  therefore,  leaves  it’s  imprint  via  the  boundary  conditions. 
The  boundary  value  problem  for  the  derivative  fields  having  been  obtained,  the  Eulerian 
derivative,  dj ,  of  the  objective  function  is  calculated  next. 
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5  The  Gradient 


Let  F(Q)  be  a  real  functional  of  the  domain,  Then  from  Eq.  (8),  its  domain  differential, 
dF( ft),  is: 

dF(ft)  =  lim  i[Fr(ftT)  -  F(ft)].  (43) 

r^O  T 

If  the  mapping:  V  dF_(tt)  is  linear  and  continuous  (in  V),  then  dF(Q)  has  a  gradient, 
Gp(fl,  V).  In  this  case,  dF( ft)  can  be  written  as: 

dF(ft)  =<  GF,V  >Y,  (44) 

where  the  duality  relation  <  •  >  is  over  Y,  the  product  space  of  Gp  and  V.  Gp  is  a 
distribution  which  is  defined  on  the  boundary,  and  can  be  considered  to  be  the  gradient 
of  F(ft).  The  objective  function  was  defined  in  Eq.  (5)  and  its  Eulerian  derivative  in  Eqs. 
(6)  and  (7).  Equation  (6)  can  be  written  as: 

r  M 

dj{ 7)  =  — Re  /  ■t/t+/Y''(A 5^>+)*5{x  —  xm)  dx.  (45) 

m= 1 

We  need  to  eliminate  the  derivative  held  in  (45).  Towards  that  let  us  define  an  auxiliary 
held,  p+,  that  satishes  the  Helmholtz  equation  below. 

M 

C+p+  =  —  ^(A 5^+)*5{x  —  xm )  x  e  ft+.  (46) 

m=  1 

In  terms  of  the  held,  p+,  (45)  for  dj{ 7)  becomes: 


dj{ 7)  =  Re  [  iP+\C+p+)  dx.  (47) 

Jn+ 

Note  that  C+  =  A  +  k+2,  and  apply  Green’s  theorem  to  (47).  After  a  straightforward 
algebra,  and  upon  simplification,  we  obtain: 


f  i/j+,(C+p+)  dx  —  f  [p+^-^+'pV\  dr  (48) 

J  0+  J  7 

£+'0+/  =  0  was  used  in  the  above  derivation.  Substituting  for  ip+'  and  t/W  from  the 
boundary  conditions  in  Eqs.  (39)  and  (40),  respectively,  we  obtain: 


n 2+ 


ip+,(£+p+)  dx  =  J  ^-P>np+  -  Pp+^j  d7+ 


'7  L 


a 


,nP+  ~  ^  'P 


,n 


d-f, 


(49) 


with 


F  —  (a  — 
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Pn  =  (a  —  ct+)V7  •  (uV7^+)  +  (a  /A2  —  a+k  2)^+v. 

In  order  to  further  eliminate  the  derivative  terms  in  Eq.  (49),  we  define  an  auxiliary 
held,  p~ ,  in  Q~,  satisfying  the  Helmholtz  equation: 


and  impose  the  conditions 


C  p  =  0,  x  G  fl  , 


p  —  p  on  7, 
P~n  =  aP+n  on  7- 


Thus,  the  auxiliary  problem  is  given  by: 


C+p+  =  —  ^^(A siJ>+)*5(x  —  xm),  x  G  fi+,  (53) 

m=  1 

=  0,  T  G  fl_,  (54) 

p_  =  p+  on  7,  (55) 

P~n  =  aP%  011  7-  (56) 

In  view  of  the  above,  we  finally  arrive  at  the  following  expression  for  dj( 7),  namely 

dj (7)  =  -Re  J  [Pp%  -  ^P,nP+^  d7,  (57) 

Reducing  the  surface  gradient  terms  by  the  Stokes’  identity  in  Eq.  (17)  then  yields: 

dj (7)  =  -Re  f  d'y  (a  -  1  )p%^%+ 

J  7 

+  (  1  —  —  ]  V7p+  •  V7^+ - —  (a~k+2  —  a+k^2)'ijj+p+  v,  (58) 

\  a  J  '  a~ 


Rewriting  Eq.  (58)  in  the  form: 


dj  =  /  Gj  v  ds , 


the  gradient,  Gj,  is  obtained  as: 


Gj  =  —Re  ( a  —  l)p+^+  +  ^1  —  —  ^  V7p+  ■  V7-0+  —  ^/c+2  —  —k  2^j  ijj+p+  (59) 
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6  Applications 


6.1  EM  Scattering 

Thus  far,  the  derivations  were  carried  out  in  terms  of  and  a1*1,  without  specifically 
referring  to  electromagnetic  (EM)  scattering.  The  EM-scattering  is  considered  now.  In 
this  case,  0  is  the  electric  (E)  or  magnetic  (H)  field,  and  a :±  is  to  be  replaced  by  e  or  /i,  the 
dielectric  and  the  permeability  constant,  respectively,  depending  upon  the  polarization 
of  the  incident  field.  It  is  well-known  [5,  34,  38]  that  an  EM-wave  propagation  can  be 
described  in  terms  of  a  TE  (transverse  electric)  and  a  TM  (transverse  magnetic)  mode 
(also  called  perpendicular  (s)  and  parallel  (p)  polarization  in  spectroscopy,  and  VV  or 
HH  in  radar  work,  see,  e.g.,  Cloude  [6]).  Transverse  implies  transverse  to  the  plane  of 
the  incidence,  that  is,  the  plane  containing  the  direction  of  the  wave  propagation  and  the 
normal  to  the  surface.  The  two  polarizations  are  schematically  illustrated  in  Figure  4.  In 
the  TE  mode,  E  is  transverse,  and  in  the  TM  mode,  it  is  H  that  is  transverse.  Referring 
to  Figure  4,  E  =  —X3E3  in  TE,  and  H  =  —X3H3  in  TM.  The  tangential  components  of 


Figure  4:  Our  Convention  of  Vector  Orientations  for  TE  and  TM  Polarizations 

both  E  and  H  are  continuous  across  the  surface.  From  this  and  Maxwell’s  equations,  the 
boundary  conditions  for  the  TE  mode  are: 


E+  = 

Etn  = 


E3  1 
E. 


3  n' 


Similar  conditions  for  the  TM  mode  are: 


Ht 


=  Hi, 

=  e+H. 


3  n’ 
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In  TE,  =  1,  whereas  in  TM,  a±  =  The  parameter,  a  =  e_/e+  =  er,  is  the  relative 
dielectric  constant.  With  these  identifications,  Eq.  (58)  for  dj  and  (59)  for  Gj  can  be 
written  down  for  the  two  polarizations.  For  the  TE  mode: 


dJrE  —  k+2(er 
Gj,te  =  k+2(er 


ijj+p+v  ds, 


l)Re  (0+p+)  . 


(60) 

(61) 


The  corresponding  expressions  for  the  TM  mode  are: 


d^JrM  —  Re  /  ds 


Or  -  i (v7 p+  ■  v7^+) 


(62) 


Gj,tm  =  Or  -  1  +  ^1  -  ^  (V7P+  ■  V7-0+)  (63) 

We  will  next  look  at  the  limiting  cases  of  the  gradient  results  for  Neumann  and  Dirichlet 
problems. 


6.2  Neumann  Surface 

Equation  (59)  can  be  used  to  obtain  the  gradient  for  both  the  Dirichlet  and  Neumann 
surface.  For  the  Neumann  condition,  the  first  term  in  the  R.  H.  S.  in  Eq.  (59),  i.e., 
(. a  —  1  )p+n'i>\1  can  be  set  to  zero  because  of  the  Neumann  boundary  condition  itself.  Next, 
letting  the  parameter,  a  — »  oo  in  the  R.  H.  S.  in  Eq.  (59)  yields: 

GJ)N  =  — R e(V7p+  •  V7^+  -  k+2^+p+). 

The  subscript  N  in  G j  implies  Neumann.  This  result  can  be  verified  by  direct  calculation 
as  follows. 

For  the  Neumann  boundary  condition,  the  unperturbed  scattering  problems  is  given 
by: 

(£+-0+)(f)  =  0,  £gD+,  (64) 

tp~(x)  =  0,  x  G  fl-,  (65) 

■0+  =  0  on  T.  (66) 

The  perturbed  problem,  i.e.,  is  also  given  by  the  above  equations.  Following  the 
procedures  in  the  preceeding  section,  the  boundary  condition  for  the  derivative  held  is 
obtained  as: 

0^  =  k+2rip+v  +  V7  •  (r:V70)  on  T,  (67) 

and  the  auxiliary  problem  for  p+  takes  the  form: 

M 

(. c+p+ )  =  -J^(A  s^+)*6(x-xm)  £efl+, 

m=  1 

p+n  =  0  on  r. 
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By  the  same  reasonings  as  above,  dj  and  Gj  for  the  Neumann  surface  are  found  to  be: 

dJN  =  —Re  f  (V7p+  •  V7V>+  —  k+2^>+p+)v  ds,  (68) 

«/7 

GJiN  =  -Re(V7p+  •  V7^+  -  k+2^+p+).  (69) 

Equations  (68)  and  (69)  are  similar  to  those  of  Norton  [36]  for  a  compactly  supported 
sound-hard  obstacle  by  variational  formulation  and  an  augmented  objective  function. 
However,  there  is  a  sign  change  in  our  results  from  Norton’s.  It  is  because  the  source 
term  in  our  auxiliary  problems  for  p+  is  opposite  to  Norton’s. 


6.3  Dirichlet  Surface 

The  same  reasonings  can  be  applied  to  the  Dirichlet  surface  also.  In  the  Dirichlet  case, 
we  neglect  the  terms  containing  and  in  Eq.  (58)  since  by  the  Dirichlet  condition, 
ijj+  is  identically  zero  on  the  surface,  Then  letting  a  — >  0  in  the  remaining  term  yields 
the  gradient  for  the  Dirichlet  surface,  namely,  G j,d  =  —Re  (V’i nP% )  •  Again,  this  result 
is  also  verified  by  direct  calculation,  as  was  done  for  the  Neumann  problem,  but  with 
,0+  =  0,  '0+/  =  —v^}\  in  the  direct  problem,  and  p+  =  0  in  the  auxiliary  problem. 
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7  Discussion 


Our  results  for  a  locally  perturbed,  infinite  surface  are  found  to  be  similar  to  the  results 
reported  in  the  literature  for  compactly  supported  objects.  This  is  not  surprising  since  the 
same  basic  equations  apply  to  either  case.  The  auxiliary  held  is  a  purely  mathematical 
construction,  and  may  not  have  a  direct  physical  meaning.  As  was  pointed  out  in  the 
text,  the  source  term  for  the  auxiliary  held,  p+,  consists  of  the  data  misfit  functions 
backpropagated  to  the  scatterer.  The  auxiliary  fields  are  then  the  fields  generated  by 
this  source  for  the  boundary  conditions  of  the  problem.  A  relation  between  the  gradient 
calculations  and  reciprocity  of  the  derivative  and  auxiliary  fields  can  thus  be  expected. 
To  illustrate  that,  let  us  rearrange  the  equations  in  a  slightly  different  form. 

Let  L  be  a  2  x  2  diagonal  matrix  whose  diagonal  elements  are  L+  and  L~ ,  and  let 
th  =  [tp+ ,  'i/j~]T ,  \E,/  =  [,0+',  ,0_/]T,  and  let  P  =  [p+,p~]T-  Moreover,  let  [q]r  denote  the 
jump  of  the  quantity,  g,  across  the  interface.  In  these  notations,  Eqs.  (l)-(4)  become: 

L\I>  =  0,  X  e  R2\t, 


Mr  =  o, 

[cy0?n]p  0. 

Similarly,  Eqs.  (37)-(42)  take  the  form: 

L\E,/  =  0,  xeR2\  r, 

W]r  =  P, 

[<#'„]  r  =  P W  > 

and  the  adjoint  held  Eqs.  (53)  through  (56)  are  similarly  recast  as: 

M 

LP  =  -Re^A  s*(xm)5(f-xm)  xeR2\T, 


m= 1 


[p]r  =  0, 

[r^p,n]p  0. 

In  these  notations,  Eq.  (58)  for  dj  takes  the  form: 

dj  =  Re  f  J  ■  n  (IT, 


(70) 


where 


J  =  —  p+[aV^']  -  W]Vp+. 


a~ 


Equation  (70)  is  the  Lorentz  reciprocity  relation  between  p+,  and  ?//  [34]. 
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Note  that  J  does  not  vanish  indicating  that  the  helds  are  not  truly  adjoint  to  each 
other.  The  non-vanishing  of  J  is  due  to  the  fact  that  the  forward  and  reciprocal  incident 
helds  do  not  act  on  the  same  scatterer,  but  the  former  on  the  unperturbed  and  the  other 
on  the  perturbed  scatterer.  In  other  words,  J  ^  0  indicates  the  variation  of  the  reciprocal 
holds.  We  have,  therefore,  used  the  word  auxiliary  instead  of  adjoint.  It  may  be  interest¬ 
ing  to  point  out  that  the  usual  picture  of  reciprocity  is  Helmholtz’s  reciprocity  in  which 
the  reciprocity  is  between  a  point  source  and  a  point  observer.  However,  Helmholtz’s 
reciprocity  does  not  always  hold  for  sources  and  observers  of  finite  extent  [21].  Lorentz’s 
reciprocity,  on  the  other  hand,  is  an  integral  relation  over  the  sources. 
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8  Conclusion 


The  inverse  scattering  solution  for  the  shape  and/or  material  parameters  of  a  scattering 
object  is  often  posed  as  a  problem  in  nonlinear  optimization  in  which  an  objective  function 
is  minimized  with  respect  to  N  (usually  large)  number  of  unknown  parameters  character¬ 
izing  the  scatterer.  The  minimization  is  often  performed  iteratively  where  the  gradient  of 
the  objective  function  in  the  parameter  vector  is  required  in  every  stage  of  iteration.  The 
application  of  the  domain  differentiation  technique  provides  a  highly  efficient  solution  to 
the  problem.  This  report  gives  a  detailed  account  of  the  domain  derivative  calculation  of 
the  gradient  for  a  one-dimensional,  locally  perturbed  infinitely  long  dielectric  surface.  Un¬ 
like  the  so-called  adjoint  field  method  which  applies  variational  procedures  to  an  objective 
function  augmented  by  the  reduced  wave  equation  via  a  Lagrange  multiplier,  the  method 
in  this  report  is  purely  algebraic  in  nature  in  that  it  evaluates  the  derivative  by  applying 
domain  differentiaion  directly  to  the  scattering  equations  without  any  other  intervening 
mathematical  formalism.  The  calculations  presented  here  are  straightforward  and  easy 
to  follow.  Moreover,  the  mathematical  transformation  of  the  scattering  problem  into  the 
corresponding  problem  for  the  differentiated  fields  can  be  visualized  explicitly.  The  so¬ 
lution  requires  an  auxiliary  scattering  problem.  The  formulation  of  and  the  motivation 
behind  the  auxiliary  problem  are  explicitly  demonstrated.  The  gradients  are  calculated 
for  the  electromagnetic  TE  and  TM  scattering  from  a  dielectric  surface  as  also  for  scale 
wave  scattering  from  a  Neumann  and  a  Dirichlet  surface.  Finally,  the  relationship  of  the 
results  thus  obtained  with  the  Lorentz  reciprocity  is  discussed. 
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List  of  Acronyms,  Abbreviations,  and  Symbols 


Acronym 

/ 

r 

k~ 

dF 

TE 

TM 


Description 

rough  interface 

Wave  number  of  upper  half  space 
Wave  number  of  lower  half  space 
Domain  differential  of  F 
Tranverse  Electric 
Transverse  Magnetic 
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